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Abstract. - Linear diamagnetism is predicted in the vortex-liquid phase of layered super- 
conductors at temperatures just below the mean-field phase transition on the basis of a high- 
temperature analysis of the corresponding frustrated XY model. The diamagnetic susceptibil- 
ity, and the Nernst signal by implication, is found to vanish with temperature as (T c o — T) s 
in the vicinity of the meanfield transition at T c o. Quantitative agreement with recent exper- 
imental observations of a diamagnetic signal in the vortex-liquid phase of high-temperature 
superconductors is obtained. 



Introduction. - The Abrikosov vortex lattice melts into an extended vortex-liquid phase 
in high-temperature superconductors subject to an external magnetic field oriented perpen- 
dicular to the conducting copper-oxygen planes that make them uppQP]- The large size in 
temperature and magnetic field of the vortex-liquid phase can be attributed to such layer 
anisotropy^S 4 5. A cross-over from a vortex-line liquid at temperatures just above the 
melting point of the Abrikosov vortex lattice to a decoupled vortex liquid at higher temper- 
ature that shows negligible correlations of the superconducting order parameter across layers 
is predicted if the vortex lattice in isolated layers melts through a continuous or a weakly 
first-order phase transition 6 . Such dimensional cross-over is observed experimentally in elec- 
tronic transport studies of the vortex-liquid phase in moderately anisotropic high-temperature 
superconductors[?|. The Abrikosov vortex lattice is predicted to sublimate directly into a de- 
coupled vortex liquid at large enough layer anisotropy, on the other hand, if the vortex lattice 
in isolated layers melts through a first-order phase transition [Sj. Electronic transport studies 
of the mixed phase in extremely layered high-temperature superconductors are consistent with 
the last sublimation scenario 

An anomalous Nernst effect is also observed in the vortex-liquid phase of high-temperature 
superconductors^. In particular, a gradient in temperature along the copper-oxygen planes 
generates an electric field perpendicular to it along the copper-oxygen planes as well. The 
low-temperature onset of the anomalous Nernst signal coincides with the melting point of 
the Abrikosov vortex lattice, while the high-temperature onset can lie above the critical tem- 
perature of the superconducting state at zero field. The authors of ref. |§] argue that this 
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effect is principally due to vortex excitations in the mixed phase of high-temperature super- 
conductors. It is then tempting to identify the cross-over between three-dimensional (3D) and 
two-dimensional (2D) vortex-liquid behavior that is predicted for layered superconductors in 
certain instances jS] with the peak in the Nernst signal. The fact that anomalous Nernst signals 
are also observed in the vortex-liquid phase of extremely layered high-temperature supercon- 
ductors that do not show the former dimensional cross-over 00 rules out that interpretation, 
however. 

The anomalous Nernst effect observed in the vortex-liquid phase of high-temperature su- 
perconductors may instead be principally due to vortex excitations in copper-oxygen planes 
that are virtually isolated from one another [3]. In this Letter, the theoretical consequences 
of that proposal are examined through a duality analysis of the uniformly frustrated XY 
model for the mixed phase of extremely type-II superconductors[S][S|. We find first that weak 
collective pinning of the vortex lattice results in a melting/decoupling temperature that does 
not extrapolate to the mean-field transition in zero field. Instead, a relatively big region of 
vortex liquid that is stabilized by random pinning centers is predicted to exist at tempera- 
tures below the mean-field transition. Second, a high-temperature expansion of the uniformly 
frustrated XY model yields linear diamagnetism at temperatures just below the mean-field 
transition. The temperature dependence of the predicted equilibrium magnetization is found 
to agree quantitatively with recent experimental reports of a diamagnetic signal extracted 
from the vortex- liquid phase of high-temperature superconductors 10 . Last, we emphasize 
that an anomalous Nernst effect is generally expected inside of the vortex liquid phased, 
where it tracks the temperature dependence shown by the diamagnetism in the vicinity of the 
mean-field phase transition. 

Vortex- Lattice Melting/Decoupling. The XY model with uniform frustration is the 
minimum theoretical description of vortex matter in extremely type-II superconductors. Both 
fluctuations of the magnetic induction and of the magnitude of the superconducting order 
parameter are neglected within this approximation. The model hence is valid deep inside 
the interior of the mixed phase. Its thermodynamics is determined by the superfluid kinetic 
energy 



which is a functional of the phase of the superconducting order parameter, e 1 ^ , over the cubic 
lattice, r. Here, J x and J y denote the local phase rigidities over nearest-neighbor links within 
layers. These are equal and constant, except over links in the vicinity of a pinning center. The 
Josephson coupling across adjacent layers, J z , shall be assumed to be constant and weak. It 
can be parameterized by J z = Jo/7' 2 , where Jo is the Gaussian stiffness of the XY model for 
each layer in isolation, and where 7' is the model anisotropy parameter. The vector potential 
Afj_ = (0, 2tt fx/a, 0) represents the magnetic induction oriented perpendicular to the layers, 
B± = $o//a 2 - Here a denotes the square lattice constant, which is of order the coherence 
length of the Cooper pairs, $0 denotes the flux quantum, and / denotes the concentration of 
vortices per site. 

The thermal/bulk average of the Josephson coupling between adjacent layers is given by 
the expression jS] 




(1) 



(cos0,,,+i> = y ]T Q(0, 1) • qr +1 (0, l^W-^-W] 



(2) 



1 



in the decoupled vortex liquid to lowest order in the fugacity yo = J z /2ksT. Here <f>ij,+\{r) = 
4>(f, I + 1) — (f)(r, I) — A z (r) is the gauge-invariant phase difference across adjacent layers I and 
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Fig. 1 - Schematic profile of the density of pinned vortices versus the total density of vortices within 
an isolated layer. Pinning centers are assumed not to crowd together. 

I + 1, and C;(l,2) = ( e -^( 1 ) e ^( 2 )) is the autocorrelation function of the superconducting 
order parameter within layer I in isolation (J z — 0). Short-range correlations on the scale of 
&d following C;(l,2) = 5oe -''i^/«2 De -#o(i) e #o(2) yields the resu lt[T2 



of the quenched inter-layer phase difference, <f>\ r = <fio(r,l + 1) — (f>o(r,l) — A z {r). It is 
set by the density of dislocations quenched into the 2D vortex lattices found in each layer at 
zero temperature in the present case of uncorrelated pinning centers. Also, above we have 
— £,2d/2 and the Josephson penetration depth A = j'a. 

In the absence of inter-layer coupling, arbitrarily weak random point pins result in a stack 
of 2D vortex lattices with dislocations quenched in^3]. Let us assume that each 2D vortex 
lattice is in a hexatic vortex glass state|14|. such that dislocations do not arrange themselves 
into grain boundaries. The quenched disorder scale 1$ that renormalizes down the interlayer 
Josephson coupling is then set by the density of such dislocations^^]- Recent theoretical 
calculations find that each isolated layer shows a net superfluid density near zero temperature 
in the collective pinning regime, where the number of dislocations quenched into each 2D 
vortex lattice is small in comparison to the number of pinned vortices |14|. Application of 
collective pinning theory to the 2D vortex lattices found in isolated layers yields a density of 
quenched-in dislocations identical to the density of Larkin domains |15| R~ 2 ~ n p (f v /vob) 2 , 
where n p denotes the density of pinned vortices per layer, where / p denotes the maximum 
pinning force, and where z/ denotes the shear modulus of the 2D vortex lattice. Here the 
critical state is assumed to be limited by plastic creep of Larkin domains by an elementary 
burgers vector b of the 2D vortex lattice. Consider now the limit of weak pinning centers that 




exp^J^l)] • exp[-#g ) +1 (2)] = e^*' 1 *. (4) 
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do not crowd together: f p — > and 7rr 2 -no <C 1, respectively, where no denotes the density 
of pinning centers per layer, and where r p denotes the range of each pinning center. Simple 
probabilistic considerations then yield the identity n p /no = 7rrp/a 2 x between the fraction of 
occupied pinning centers and the ratio of the effective area of each pinning center to the area 
per vortex, a 2 x = a 2 //. This yields the result n p — (no • 7rr 2 )n vx for the density of pinned 
vortices|16|. where n vx = l/ a vx i s the density of vortices per layer. Finally, substitution of the 
estimate Vq = (7r/4)n vx J for the shear modulus [T7) yields the result R~ 2 ~ (f p r p / J ) 2 n for 
the density of Larkin domains, which is independent of magnetic field. Note, however, that 
all of the above is valid only in the 2D collective pinning regime that exists at perpendicular 
magnetic fields above the threshold Bc P D ^ ~ (/ p / Jo) 2 $o, in which case many vortices are 
pinned in each layer within a Larkin domain of dimensions R c x R c |14|. Single- vortex pinning 
exists at magnetic field below that threshold, on the other hand, in which case each Larkin 
domain contains only a single pinned vortex: n p = R~ 2 . Assembling the above suggests the 
profile for the density of pinned vortices per layer versus the density of vortices that is depicted 
by n g- HI K implies that the quenched disorder scale 1$ ~ R c is independent of magnetic fields 

above the threshold B^ ^ ~ <£> /-R 2 - 

We are finally in a position to determine the melting/decoupling line of the 3D vortex lat- 
tice at temperatures outside of the 2D critical regime, ^2D ~ a vx , at big enough perpendicular 
magnetic fields such that Larkin domains can be defined, B± 3> B^ D ^ . The identification of 
the separation between dislocations quenched into each 2D vortex lattice with the 2D Larkin 
scale I,), ~ R c , necessarily yields the inequality <C 1$ in such case. At temperatures 
lying inside of the interval [Tm D \ J 1 ] 215 -*] bounded by melting of the 2D vortex-lattice and by 
the Kosterlitz-Thouless transition in isolated layers, yet lying outside of the 2D critical regime, 
a partial duality analysis of the pristine layered XY model with uniform frustration Q finds 
a first-order melting/decoupling transition of the 3D vortex lattice at interlayer Josephson 
coupling 5 

(cos0 M+1 ) ~ 1/2. (5) 

The first-order nature of this melting/decoupling line and its coincidence with the contour 
defined above is consistent both with Monte Carlo simulations of the same modeljTT] and 
with elastic medium descriptions of the vortex lattice in layered superconductors 4 . Observe 
now that the criterion (JjjJ for first-order melting/decoupling should remain valid in the present 
regime of weak pinning such that 3> Substitution of expression for the inter-layer 
"cosine" in the decoupled vortex liquid then yields the melting/decoupling field 

Bn={^-\[^? ] )\ (6) 

where B^} ~ g 2 (Jo/fc,BT)(<I>o/A[)) is the melting/decoupling field in the pristine limit |4*]|5)|11|. 
These results are summarized by the phase diagram shown in fig. [21 The short sections 
of dashed and solid lines that emanate perpendicularly from the horizontal axis originate 
respectively from the decoupling cross-over © and the second-order phase transition shown 
by the layered XY model in the absence of uniform frustration (cf. ref. ^H|)- We conclude this 
section by observing that the melting/decoupling line does not extrapolate to the mean-field 
critical temperature at zero-field [Jo(T c q) = 0] due to the presence of dislocations quenched 
into the weakly pinned vortex lattices found in isolated layers. 

Vortex- Liquid Diamagnetism. The phase diagram for the mixed phase of layered 
superconductors shown by fig. implies a large region of vortex liquid in the vicinity of 
the meanfield transition at zero magnetic field because of the effects of random point pins. 
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Fig. 2 - Schematic phase diagram near the critical temperature at zero magnetic field. Continuity with 
the zero-field limit of the XY model Q implies a lower critical point, which is observed experimentally 
in the phase diagram of high-temperature superconductors with point disorder. (See ref. 

In particular, the equilibrium diamagnetic susceptibility due to the emergence of Cooper 
pairs is well defined at temperatures inside of the window [T c , T c q]. The former quantity- 
can be obtained from the uniformly frustrated XY model JQ) in the vicinity of the mean- 
field transition via a high-temperature expansion in powers of the fugacity zo = J/2ksT 
In particular, a duality analysis yields that the corresponding partition function is 
approximated by Zxy — Zq + Z4 as zq — > 0, where Zq = Yluj) Io(Jij I^bT), and where 
Z^/Zq — 2cos(27r/)^pJJ/ i .v g pti(J<j/AsBT). Here t\(x) = 1% (x)/Io (x ) is the ratio between 
a first-order and a zero-order modified Bessel function, which is approximately ti(x) = x/2 
for \x\ <C 1. Also, (ij) represents nearest-neighbor links within layers, and □ represents 
elementary plaquettes within layers. The equilibrium magnetization is given by Mj_ = 
—d(GxY/V)/dB^ in the extreme type-II limit, where Gxy — ~k B T ■ \nZ X Y is the Gibbs 
free energy, and where V is the volume. Substitution of the previous high-temperature ap- 
proximation yields 

M ± = -(2V$o)(J4/d)(J4/2fc s T) 3 sin(27r/), (7) 

where J4 = (Il(ij)eP ^w) ly/4 > an d where d denotes the spacing between layers [27)]. The mag- 
netization therefore varies linearly with vanishing magnetic field like M± = xH±, with a 
diamagnetic susceptibility 

An X - ~ K - 2 (Jo/2k B T) 3 (J 4 /J ) 4 (a/0 2 - (8) 

Here k = is the usual ratio of the London penetration depth to the coherence length 

of the Cooper pairs. The former is related to the Gaussian phase stiffness of each layer by 
Jo - $oW16^ 3 A|. 

Non-linear diamagnetism is observed experimentally in the vortex liquid phase of the 
extremely layered high-temperature superconductor L^S^CaC^Os+a (BSCCO), at temper- 
atures just above the superconducting transition in zero field |2JJ • Linear diamagnetism is 
displayed at yet higher temperature in the same samples, on the other hand. By Eq. JJJ, 
the uniformly frustrated XY model JJJ predicts such linear diamagnetism, M± = x^±, at 
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Fig. 3 - Fit of diamagnetic signal extracted from underdoped BSCCO (ref. |10| ) to the theoretical 
prediction near the mean-field transition [Eq. @]. 

perpendicular magnetic fields that are small compared to the upper-critical scale <I>o/27ra 2 , 
at temperatures just below the mean-field phase transition. The corresponding diamagnetic 
susceptibility predicted by the XY model is given by Eq. JSJ. Use of the relation quoted 
previously between the Gaussian phase rigidity within planes and the London penetration 
length in conjunction with physical parameters n = 100, A^(0) = 0.2 /im, and d = 1.5 nm 
appropriate for BSCCO [22] yields the estimate 



for the diamagnetic susceptibility of that material in the vicinity of the mean-field phase 
transition. Here n s (T)/n = )? L (Q) / \ 2 L (T) is the superfluid fraction. The mean-field superfluid 
density expected from a pristine <i-wave state in 2D is approximately 2/3 the corresponding s- 
wave result in the vicinity of the meanfield transition at zero field 23]; i.e., n s /n = (4/3)(l — t), 
where t — T/T c o is the reduced temperature. Equation © then implies that the diamagnetic 
susceptibility vanishes like (1 — t) 3 with temperature as it approaches the mean-field transition. 
Figure [3] displays the cube-root of the diamagnetic signal extracted experimentally in ref. |10| 
from an underdoped sample of BSCCO with T c = 50° K, in perpendicular magnetic field 
H± = 14 T, as a function of temperature. The solid line is a fit to the linear diamagnetism, 
M± = X-H-Li predicted by the high-temperature expansion of the uniformly frustrated XY 
model, Eq. @, with n s /n = (4/3)i(l - t), T c0 = 158° K, and with XY model parameter 
a = 0.28 ( Jo/ Ji) 2 £,- The success of the fit indicates that the onset of the diamagnetic signal 
observed in the vortex liquid phase of high-temperature superconductors reflects nothing other 
than the mean-field phase transition at which Cooper pairs emerge. The large suppression of 
T c compared to the meanfield transition temperature T c o obtained here can be accounted for 
by quenched disordering of the superconducting order parameter, which could be generic to 
under-doped high-temperature superconductors 

Anomalous Nernst Effect and Conclusions. - A gradient in temperature along the layers 
in the vortex liquid phase of high-temperature superconductors generates a voltage in the 
perpendicular direction within the layers^-. In particular, the Nernst signal defined by the 



4ttx/mo 



[(503 o K/T)(n s /n)] 3 (J 4 /J ) 4 (a/O 2 A/Tm 



(9) 
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ratio e y — E y /d x T between the electric field that is generated and the gradient in temperature 
peaks inside of the vortex liquid. Standard transport theory yields the identity|U] 

e y =p x - a s xy (10) 

between the Nernst signal and the product of the flux-flow electrical resistivity p x with the 
off-diagonal Peltier coefficient a s xy . Also, application of Ginzburg-Landau theory for the su- 
perconducting order parameter yields the estimate 25 

<y = 0M ± (11) 

for the Peltier coefficient near the mean- field transition, where (3 is of order T . Observe 
now that the flux-flow resistance increases with temperature in the vortex liquid, while the 
equilibrium magnetization decreases with temperature there [cf. Eq. JSJ]. Substitution of 
the estimate 111 If into the identity IjlUfl then yields (i) that the low-temperature onset of 
the anomalous Nernst signal is given by the melting/decoupling temperature of the vortex 
lattice. Also, the linear diamagnetism © extracted from the high-temperature regime of the 
frustrated XY model implies (ii) that the anomalous Nernst signal vanishes with temperature 
at the mean- field transition as (T c q — T) 3 . Where exactly the Nernst signal peaks inside of 
the vortex- liquid phase depends on how pinning affects the flux- flow resistance [22J, which is 
beyond the scope of the paper. 

In conclusion, a high-tcmperature analysis of the layered XY model with uniform frus- 
tration finds that the simultaneous onset of linear diamagnetism and of an anomalous Nernst 
effect in the normal phase of high-temperature superconductors ^2] can be identified with the 
mean-field transition for Cooper pairing. The low-temperature onset of the anomalous Nernst 
signal at the melting/decoupling line of the vortex lattice was also found to be depressed sub- 
stantially by the presence of dislocations quenched into the vortex lattice in isolated layers. 

* * * 
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